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o : 

' The extended Gari-Kriimpelmann (GK) model of nucleon electromagnetic form fac- 
tors, in which the p, p', w, w' and vector meson pole contributions evolve at high 
momentum transfer to conform to the predictions of perturbative QCD (pQCD), 

\ was recently shown to provide a very good overall fit to all the nucleon electromag- 

Tj" ' netic form factor (emff ) data, including the preliminary Rp and Rn polarization data 

' available in 2002, but excluding the older Gep and Geu differential cross section data 

^SJ ■ that was inconsistent with the Rp and Rn data. The recently published final version 

. of the polarization data of the electric to magnetic ratio Rp differs little from the 

\ preliminary values for the former, but the new values of Rn are midway between 

I the preliminary values and those inferred from the differential cross section data and 

f-*) ' the Rosenbluth separation. A new fit of the parameters of the same model has been 

■ made with the final Rp and Rn data replacing the preliminary values and the addi- 

. tion of some new Rn and Gmu data. Again there is a good fit to all the data when 

\ excluding the differential cross section data that is inconsistent with the polarization 

H \ data. This includes a very good fit of the Rp data, which was not possible when the 

^ ' differential cross section G(^En) data was used in place of the polarization Rn data. 

C ■ Thus the change between the preliminary and final Rn data, while substantial, has 

^ . not impeded the good simultaneous fit to the neutron and proton data. The param- 

^ , eters, fit to the data and predictions of the new model are compared to those of the 

\ previous models. Low momentum structures that appear in some data analyses are 
partially reproduced by the model. 

PACS numbers: 13.40.Gp, 21.10.Ft 



I. INTRODUCTION 

The Gari-Kriimpelmann (GK) model of nucleon electromagnetic form factors 0| incor- 
porates the physical features of meson exchange (in the vector dominance approximation) 
at low momentum transfers, and asymptotic freedom at high momentum transfers on the 
QCD scale. In previous papers and j3| the GK model was extended in its parameteriza- 
tion of the hadronic form factors, by including the width of the p meson, and adding p' and 
uo' exchange. When the polarization data for the ratios -R„ and Rn became available it was 
shown that the GKex(02) model (Eqs. 5, 6 and 7 of 13) provided a particularly good 
fit to the neutron and proton electromagnetic form factor data when polarization values 
for Gep and Geu were used, replacing the inconsistent values obtained from the differential 
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cross sections. 

It is shown here that the quahty of the fit is maintained, and the fitted parameters are 
closer to the physical values, when the revised final vesions of the Rp and Rn data, and some 
new data, are used. It is also noted that relatively narrow structures appearing in some data 
near ^ Q.2{GeV / (?) are partially reproduced by the data. The model parameters were 
negligibly influenced by the small weight of the data in the region of the structures. Instead 
the effect can be traced to the influence of the full momentum range of the data on the sign 
variation of the residues of the vector meson poles. 



II. THE NUCLEON EMFF MODEL 

In fitting the nucleon emff data including the newly available Rn 0] and Rp jH results we 
use the extended GK model GK(02) of Ref. as defined by Eqs. 5, 6 and 7 of 's'], which 
has the following features: 

(1) The photon- nucleon interaction is mediated by the exchange of the vector- mesons p, 
cj, 0, p'(1.45 GeV) and cj'(1.419 GeV). 

(2) The vector mesons are treated in the zero-width (pole) approximation, with the ex- 
ception of the p meson. 

(3) An approximation to the dispersion integral, that incorporates the width, is used for 
the p meson. 

(4) It uses the QCD cut-off A2 for the helicity flip meson-nucleon form factors, rather than 
the meson-nucleon cut-off Ai used in previous versions. 

(5) The evolution of the logarithmic dependence on is controlled by the quark-nucleon 
cut-off Ad, along with Aqcd- The latter is fixed in its experimental range at 0.15 GeV. 

(6) Asymptotically the form factors have the pQCD behavior (see Eq. 8 of 3]). 

(7) The form factor Fi{Q'^) vanishes at = 0, and it and also Ffi^Q"^) decreases more 
rapidly at large than the other meson form factors. This conforms to the Zweig 
rule imposed by the ss structure of the meson 



III. DATA BASE AND FITTING PROCEDURE 

In the fit GKex(02S) of |3| the data of Refs.[7-14, 16-36] in that paper was used, with 
the omission of the Gep values for > 1.75 GeV^/c^ of that Ref. [7] and the Geu values 
for Q2 > 0.779 GeVVc^ of Refs.[9, 17, 18] there. In the present fit we substitute the Rp 
values of and the Rn values of Q for the preliminary values of the Ref[5] and Ref[14] of 
the previous paper As noted above, the Rp values change little (although there is one 
more data point), while there is a substantial change in the Rn values. In addition we add 



and Gmu data of 0- 
3| the 8 photon coupling constants of the four vector mesons 
treated as poles and 4 scale parameters were varied. However the present fit was made with 
fixed N=l (implying negligible error in the p meson dispersion relation evaluation) and with 



the newly available Rn data of 
As in the case GKex(02S) of 
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Aqcd fixed at the physical value of 0.15 GeV/c. The 12 free parameters were optimized 
using a Mathematica program that incorporates the Levenberg-Marquardt method. 

There are two pairs of parameters which are so strongly correlated that the fit varies 
negligibly over large ranges of the parameters when both vary in a correlated manner. In 
fits to the data the photon-p' coupling, gpij jp'^ is small. The contribution of the p' to F\ is 
negligible, but if Hpi is large the p' contribution to F2 is of importance and best fits the data 
when the product of the two parameters is near .08 . For Kpi running from 12.0 to infinity, 
Qp'l fp' runs from .007 to zero. For this whole range the value only changes by < 0.1. We 
constrain Kpi to the smallest value in this range with the correlated largest value of Qp'/ fpi- 

The other strongly correlated parameter pair is /t^ and p^ which appear together in 
the hadronic form factor F2 . In the fits /i^ is small (< Q.2GeV / c). Consequently for 
Q> Q.hGeV/c, F2 is approximately proportional to k^/ p^. As /x^ runs from 0.2 to zero, 
runs from 0.01 to zero and the best fit again only changes < 0.1. Constraining both p^ to 
0.2 and k.^), to 0.01 only increases the best fit by 0.06. Only the remaining 8 parameters 
are individually sensitive to the data fit. 

The resulting parameters and the quality of the fit to this data set, GKex(05), is compared 
with that of the GKex(02S) and GKex(Ol) models in Tables I and II. The fit of these models 
to the complete data set (including the differential cross section data for electric form factors 
inconsistent with polarization data) is displayed in Figs. 1-6. The model GKex(Ol), which 
did not include the uj' meson, was fitted to all the differential cross section data, including 
the available Rp data but in the absence of the data. The present GKex(05) model differs 
from GKex(02S) only due to the substitution of the newer polarization data and inclusion 
of the few new Rn and Gun points. 

IV. RESULTS 

Table H] presents the parameters which minimize for the above 3 cases, to their own 
data sets For all 3 parameter sets the hadronic form factor cut-off masses, Ai, A2, A15, and 
P(l). are reasonable. The relatively large value of A2, which controls the spin-flip suppression 
in QCD, is consistent with the slow approach to asymptopia observed in polarized hadron 
scattering. For the case in which Aqcd is a fltted parameter, as well as the two for which it 
is flxed, it is consistent with high energy experiment. The addition of the ci;'(1.419) meson in 
GKex(02S) and especially in GKex(05) has moved closer to the expected small negative 
value than all earlier flts, but the remaining small positive value implies some contribution 
may be required from a still higher mass isoscalar meson. However the adequacy of the flts is 
an indication that form factors with more poles would produce similar flts to those already 
obtained. 

In Table im the values of for the present data set are listed for the three cases, detailing 
the contribution from each of the six form factors. The values for Gmp, Gep, Rp and Rn 
are less than the number of data points, but for Gmu, and to a lesser extent for Geu, the 
values are substantially larger than the number of data points. These large values of x^ are 
caused by the inconsistent fluctuations of the data from point to point for <2 GeV^/c^. 
Fitting the data better in that range would require unreasonable oscillations in a model. 

Fig. 1 shows the substantial effect on the best flt model of the new highest momentum 
transfer polarization i?„ value 01 • The experimental values and the model predictions are 
substantially larger than the phenomenological Galster prediction, but far less than the 
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FIG. I: Rn, the ratio ^J-nGEn/GMn- Comparison of the models GKex(Ol) [solid], GKex(02SHdash- 
dotted] and GKex(05) [dotted] with the data. The dashed cm^ve is R^"^'^'"'' (Q^) of Eq. 1 of Q- The 
experimental points are described in the text |^,|^. The data symbols are listed beside the figure. 



values determined from the differential cross section data. The validity of the model up to 
3.5 GeV'^/c^ will be tested when recent JLab data [8] has been analyzed. However, as seen 
in Fig. 2, the substantial change in the model parameters required to fit the new Rn data 
(together with new Rp and Gmu data) has only a small effect on the fit to the Rp data. This 
is remarkable when compared to the much poorer fit to the Rp data required by fitting the 
GKex(Ol) data set, which differs from the Rn GKex(Ol) curve by about the same amount 
as the latter differs from the GKex(02S) curve. As shown in Figs. 3-6 there is little change 
in the Gep and Gmp predictions up to large momentum transfers, while Gmu and Geu are 
substantially altered at the larger momentum transfers. 




FIG. 2: Rp, the ratio fJ.pGEp/GMp- (a) Comparison of the models GKex(Ol) [solid], GKex(02S) 
[dash-dotted] and GKex(05) [dotted] with the data. The polarization data is from [51, and from 
Ref. [12] of (b) Comparison of model GKex(05) [dotted] with nth new data of [ll]. The data 
symbols are listed beside the figure. 



It is noteworthy that, while the parameters of this model were fitted to the whole mo- 
mentum transfer region of the available data, the model reproduces, at least in part, narrow 
structure observed at very low momentum transfer, between 0.1 and 0.4 GeV'^/c^. This was 
noted in the data available in 2003 by , and may be confirmed by recent data taken with 
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FIG. 3: i(a) Geti normalized to Gd- Comparison of the models GKex(Ol) [solid], GKex(02S) [dash- 
dotted] and GKex(05) [dotted] with the data of Refs. [9, 17, 18. 26, k 28-36] of B. The data symbols 
are listed beside the figure. The points labelled by open circles are obtained by multiplying Rn 
data 01 and by the Gmu of GKex(05) normalized by /J^nGd (b) Expansion of the range to 
Q2 < 2.0 GeV^/cP. 
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FIG. 4: Gep normalized to Gd- Comparison of the models GKex(Ol) [solid], GKex(02S) [dash- 
dotted] and GKex(05) [dotted] with the data of Refs. [6-11] of 0|. The data symbols are listed 
beside the figure. The points labelled by open circles are obtained by multiplying Rp data of |^ 
and Ref. 13 of 0| by the Gmp of GKex(05) normalized by fipGd (a) The full data range, (b) 
Expansion of the range < 2.0 GeV^/c^. 



the BLAST detector at Bates Laboratory [9^. In [10| the structure is phenomenologically 
related to a pion cloud. In a VDM model, such as the ones discussed here, the pion cloud is 
represented by pion pairs and triplets clustered into vector mesons. Figs. 3-6 (b) show the 
data in the expanded scale of < 2.0 GeV^/c^ to better observe the small structures in the 
data and models. For both Gmp in Fig. 5(b) and Gmu in Fig. 6(b) a minimum appears near 
0.2 GeV^/c^ in the published data and is reproduced by the models of this article. For 
Gep the data also shows a minimum in the same vicinity, but our models. Fig. 4(b), only 
have a horizontal inflection there. Our models show no low momentum-transfer structures 
for Geu, Fig. 3(b), and the presence of the maximum indicated in jl^j is not clear in the 
published data. The recent experiments may provide a clarification. The proton polar- 
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(a) (b) 

FIG. 5: Gmp normalized to ^pG^- Comparison of the models GKex(Ol) [solid], GKex(02S) [dash- 
dotted] and GKex(05) [dotted] with the data of Refs [6-9, 11, 12, 19 & 20] of 3], . The data symbols 
are listed beside the figure, (a) The full data range, (b) Expansion of the range < 2.0 GeV^/c^. 
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FIG. 6: Gmti normalized to UnGd- Comparison of the models GKex(Ol) [solid], GKex(02S) [dash- 
dotted] and GKex(05) [dotted] with the data of and of Refs. [9, 16-18, & 21-27] of Q. The 
data symbols are listed beside the figure, (a) The full data range, (b) Expansion of the range 
Q2 < 2.0 GeVVc^. 



ization data has now become available [llj] , too late to be included in the model parameter 
fitting. That data is compared with model GKex(05) in Fig. 2(b), showing two data points 
(out of eight) that are lower than the model curve by about i.5 standard deviations. To the 
extent that the model does have structure near ^ 0.2 GeV^/c^ the effect can be traced 
to the opposite signs of the u and uj' meson contributions to combined with the more 
rapid decrease of the uj contribution in that region. 



V. CONCLUSIONS 



The incorporation of data that has become available since the publication of p|, espe- 
cially the Rn data of [4], was used to refit the parameters GKex(02S), resulting in the version 
GKex(05). Only the predictions for Rn and Geu at higher momentum-transfers were substan- 
tially changed, and a good fit to the data maintained. The predicted values of Rn continue 
to rise above the Galster curve with increasing momentum transfer, a prediction which may 
soon be checked by new data 0. The fitted parameters agree with known constraints and 
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the model is consistent with VMD at lower momentum-transfers, while approaching pQCD 
at high momentum-transfers. 

It also has been noted that the predictions of the model include at least some of the 
structure recently noted near ?a 0.2 GeV^/c^. 
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TABLE I: Model parameters. Common to all models are = 3.706, Kg = —0.12, rup = 0.776 GeV, 
= 0.784 GeV, = 1.019 GeV, = 1.45 GeV and m^^ = 1.419 GeV. 



Parameters 


Models 




GKex(Ol) GKex(02S) 


GK(05) 


dp'/fp' 


0.0636 


0.0401 


0.007208 


Kp' 


-0.4175 


6.8190 


12.0 


9w/fuj 


0.7918 


0.6739 


0.7021 




5.1109 


0.8762 


0.4027 


9<t>/U 


-0.3011 


-0.1676 


-0.1711 




13.4385 


7.0172 


0.01 




1.1915 


0.8544 


0.2 






0.2552 


0.164 






1.4916 


-2.973 


Ai 


0.9660 


0.9407 


0.93088 


Ad 


1.3406 


1.2111 


1.181 


A2 


2.1382 


2.7891 


2.6115 


Aqcd 


0.1163 


0.150 (a) 


0.150 (a) 


N 


1.0 (a) 


1.0 (a) 


1.0 (a) 



(a) not varied 



TABLE II: Contributions to the standard deviation, for each form factor, from the present data 
set (Section 3) for each of the three models. The number of data points contributing is in the 
second column. The polarization data is only used for the Rn,p values and only differential cross 
section data is used for the GMp.Ep,Mn,En values (and cold neutron scattering for the G£;„ slope at 
Q=0). 



Data 
type 


Data 
size 


Models 

GKex(Ol) GKex(02S) GKex(05) 


Gmp 


68 


43.3 


47.9 


51.5 


Gep 


36 


26.9 


30.5 


28.4 


Gmu 


39 


127.8 


129.4 


124.9 


Geu 


15 


25.3 


24.1 


25.1 


Rp 


22 


135.8 


10.4 


10.3 


Rn 


5 


1.9 


5.6 


1.1 


Total 


185 


361.0 


248.0 


241.3 



